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ABSTRACT 

Most theoretical models invoke quasar driven outflows to quench star formation in massive 
galaxies, and this feedback mechanism is required to account for the population of old and 
passive galaxies observed in the local universe. The discovery of massive, old and passive 
galaxies at z~2, implies that such quasar feedback onto the host galaxy must have been at work 
very early on, close to the reionization epoch. We have observed the [CII]158^m transition 
in SDSSJ1 14816.64+525150.3 that, at z=6.4189, is one of the most distant quasars known. 
We detect broad wings of the line tracing a quasar-driven massive outflow. This is the most 
distant massive outflow ever detected and is likely tracing the long sought quasar feedback, 
already at work in the early Universe. The outflow is marginally resolved on scales of ~16 kpc, 
implying that the outflow can really affect the whole galaxy, as required by quasar feedback 
models. The inferred outflow rate, M > 3500 M Q yr _1 , is the highest ever found. At this rate 
the outflow can clean the gas in the host galaxy, and therefore quench star formation, in a few 
million years. 
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1 INTRODUCTION 

The old stellar population and low gas content characterizing 
local massive galaxies, as well as their steeply declining num- 
ber at high masses, have been difficult to explain in most theo- 
retical scenarios of galaxy formation. Indeed, in the absence of 
"negative feedback" on star formation (i.e. a mechanism to de- 
crease star formation), theoretical models expect a much larger 
number of massive galaxi es than observed and characterized b y 
young stellar populations (Menci et al. 20061; iBower et all l2006h . 
To reconcile theory with observations, most models invoke mas- 
sive outflows driven by the radiation pressure generated by quasars, 
which are expected to clear their host galax ies of the bulk of 



their gas, therefore quenching star formation (Silk & Rees 


1998; 


Granato et al. 


2004: Di Matteo et alj 20051; Sprinsel et al. 


2005; 


Croton et al. 


20061: Narayanan et all 120081: iHopkins et all 


2008; 


Fabian et alj|2009l; lKinsll20ld: iKins et al]|201ll: 1 Hopkins & Elvis! 


2010t Zubovas & Kindl2012h. Evidence for such negative feed- 



back has been found only recently in quasar hosts through the 
discovery of massive outflows, in particular t hrough the detection 
of P-Cygni profil es of far-IR OH transitions (Fische ret al.ll2O10l ; 
ISturm et all 1201 ll) . the detectio n high velocity wings of molecu- 
lar emission lines (CO, HCN, iFemglio et al7l2QloflAlatalo et alj 



l201ll;lAalto et alj|2012h and the detection of high velocity neutral 
gas in absorption and of high velocity ionized ga s (Nesvad ba et alj 
l2010l : iRupke & Veilleuxll201ll ; iGreene et alj|201 lb . Evidence for 
quasar-driven massive outflows have also been found through the 
bright [CII]158yum fine structure line; indeed, Herschel spectra of 
nearby quasars have revealed prominent [CII] broad wings, asso- 
ciated with the molecul ar outflow traced by the OH absorption 
lines dFischeretal1l2010l. Sturm et al. in prep.). Quasar-driven out- 
flows have been det ected also at high redshift (although exploit- 
ing different tracers , lAlexander et aljEoiOt iNesvadba et al.ll201 it 
iHarrison et al.l2012h . Moreover, evidence has been found that these 
outflows are indeed quenching star formation in their host galaxies 
dCano-Dfaz et all2012l ; iFarrah et al]|2012l; ISteinhardt & Silverman] 
l201ll ; |Page et alj|2012l) . as expected by models. 

However, the disco very of old (age > 2 Gyr), passiv e and 
massive galaxies at z~2 (cimatti et alj | 2004l; ISaracco et al.1120051 ; 
iMcCarthv et alJl2004l;lGlazebrook et alJl2004T) . when the Universe 
had an age of only ~3 Gyr, requires that the quasar feedback 
quenching mechanism must have been at work already at z>6 
(age of the Universe less than one Gyr), i.e. close to the re- 
ioniza tion epoch. Quasar driven winds have been observed up to 
z~6 dMaiolino et alj 1200 ll . 12004). however these are associated 
with ionized gas in the vicinity of the black hole, accounting only 



© 2002 RAS 



2 Maiolino et al. 



for a tiny fraction of the total gas in the host galaxy. So far, no ob- 
servational evidence was found for the massive, quasar driven out- 
flows at z>6 required by feedback models to explain the population 
of old massive galaxies at z~2. 

Here we focus on one of the most distant quasars known, 
SDSSJ1 14816.6 4+525 150.3 (hereafter Jl 148+5251), at z=6.4189 
dFan et all [2003). CO observations have revealed a large reser- 
voir of mol ecular gas, Mh, ~ 2 x 10 10 Mp, in the quasar 
host galaxy dWalteret al]|2003l : iBertoldi et alll2003bl) . The strong 
far-IR thermal emission inferred from (sub-)millimeter observa- 
tions reveals vigorous star formation in th e host galaxy, wit h 
SFR ~ 3000 M Q yr 1 dBertoldi et al]|2003ak iBeelen et alfl2006T) . 
Jl 148+5251 is also the first high redshift galaxy in wh ich the the 
[CII]158yum line was discovered jMaiolino et"atl f2005). High res- 
olution mapping of the same line with the IRAM PdBI revealed 
that most of the emission is confined within ~ 1.5 kpc, indicating 
that m ost of the star form ation is occurring within a very compact 
region dWalter et al.ll2009t) . 

Previous [CII] observations of Jl 148+5251 did not have a 
bandwidth large enough to allow the investigation of broad wings 
tracing outflows, as in local quasars. In this letter we present new 
IRAM PdBI observations of Jl 148+5251 that, thanks to the wide 
bandwidth offered by the new correlator, have allowed us to dis- 
cover broad [CII] wings tracing a very massive and energetic out- 
flow in the host galaxy of this early quasar. We show that the prop- 
erties of this outflow are consistent with the expectations of quasar 
feedback models. 

We assume the concordance A-cosmology with Hp = 
70.3 km s -1 Mpc" 1 , Q A = 0.73 and £l m = 0.27 dKomatsu etaT] 

I201I . 



2 OBSERVATIONS AND DATA REDUCTION 

Observations with the IRAM PdBI were obtained mostly in April 
2011 in D configuration (mostly with 1.5 < PWV < 3.5 mm), 
while a few hours where also obtained in January 2011 in C+D 
configuration (PWV < 1.5 mm). The resulting synthesized beam 
is 2.2" x 1.8". The receivers were tuned to 256.172 GHz, which 
is the rest frame frequency of [C II] at the redshift of the quasar, 
z=6.4189 dMaiolino et ail [2005). The following flux calibrators 
were used: 3C454.3, MWC349, 0923+392, 1150+497, 3C273, 
3C345, 1144+542, J1208+546, J1041+525, 1055+018. Uncer- 
tainties on the absolute flux calibration are 20%. The total on- 
source integration time was 17.5 hours, resulting in a sensitivity 
of 0.08 Jy km s _1 beam 1 in a channel with width 100 km s -1 , 

The data were reduced by using the CLIC and MAPPING 
packages, within the GILDAS-IRAM software. Cleaning of the re- 
sulting maps was run by selecting the clean components on an area 
of ~ 3" around the peak of the emission. For each map the resulting 
residuals are below the lcr error, ensuring that sidelobes are prop- 
erly cleaned away. Anyhow, as discussed in the following, the size 
determination has been investigated directly on the uv data, there- 
fore independently of the cleaning. 



3 RESULTS 

3.1 Detection of broad wings 

The continuum was subtracted from the uv data by estimating 
its level from the channels at v < -1300 km s -1 and at v > 
+ 1300 km s . The inferred continuum flux is 3.7 mjy, which is 
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Figure 1. IRAM PdBI continuum-subtracted spectrum of the [CII]158//m 
line, redshifted to 256.172 GHz, in the host galaxy of the quasar 
Jl 148+5152 extracted from an aperture with a diameter of 4", top, and 
6", bottom. The spectrum has been resampled to a bin size of 85 km s . 
The red lines show a double Gaussian fit (FWHM=345 km s -1 and 
FWHM=2030 km s ) to the line profile, while the blue line shows the 
sum of the two Gaussian components. 



fully consistent with the value expected (4 mJy) from the bolomet- 
ric observations dBertoldi et al.ll2003ah . once the frequency range 
of the latter and the steep shape of the thermal spectrum are taken 
into account. 

Fig- Eh shows the continuum-subtracted spectrum, extracted 
from an aperture of 4" (corresponding to a physical size of 1 1 kpc). 
Fig- Hp shows the spectrum extracted from a larger aperture of 6" 
that, although noisier than the former spectrum, recovers residual 
flux associated with the beam wings and with any extended com- 
ponent. 

The spectrum shows a clear narrow [CII]158yum emis- 
sion line, which was already detected b y previous observations 
dMaiolinoet ai]|2005l : IWalter et alj|2009l) . However, thanks to the 
much wider bandwidth relative to previous data, our new spectrum 
reveals broad [CII] wings extending to about ±1300 km s _1 . These 
are indicative of a powerful outflow, in analogy with the broad 
wings that have been observed in the molecular and fine struc- 
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ture lines of local quasar host galaxies. The spectrum is fifl with 
a narrow (FWHM=345 km s -1 ) and a broad (FWHM=2030 km 
s -1 ) Gaussian, as shown in Fig. [T] resulting into^Lj = 1.11. By 
removing the broad component the x 1 increases from 195 to 229 
(with 175 degrees of freedom), implying that the broad component 
is required at a confidence level higher than 99.9%. 

Theoretical models predict that starburst-driv en winds can- 
not reach velocitie s higher than 600 km s" 1 dMartinl 120051 : 
Thack er et al. 2006), therefore the high velocities observed in the 
[CII] wings of Jl 148+525 1 strongly favor radiation pressure from 
the quasar nucleus as the main driving mechanism of the outflow. 
Quasar radiation pressure is favored, relative to SN-driven winds, 
also based on energetics arguments, as discussed later on. 



3.2 Outflowing gas mass 

The luminosity of the broad [CII] component allows us to infer a 
lower limit of the outflowing atomic gas mass, using, 

Moutu (atomic) = Q ?? ^ 0.7 L [C1I] ^ / 1.4 lO -4 ^ 1 + 2 e -911 ^ + Wn 



M G -■ \ Lo )\ Xc + / 2e- 9 ' K / T 

faailev-Dunsheathet ai] |2010l) . where X c + is the C + abundance 
per hydrogen atom, T is the gas temperature, n is the gas den- 
sity and n crit is the critical density of the [CII]158/im transition 
(~ 3 x 10 3 cm -3 ). By assuming a density much higher than the 
critical density, Eq.Q]gives a lower limit on the mass of atomic gas. 
The quasar-driven outflows observed locally are characterized by a 
wide range of de nsities, including both dense clumps and diffuse, 
low d ensity gas dCicone et alj2012l : lAalto et alj2012l : lFischer et alJ 
2010), therefore our assumption on the gas density gives a con- 
servative lower limit on the outflowing gas mass. We assume a 
temperature of 200 K, however the dependence on the tempera- 
ture is weak (going from 100 K to 1000 K the implied gas mass 
is within 20% of the value obtained at 200 K). Furthermore, we 
assume a C + abundance ty pical of PDRs, i.e. X c + = 1.4 x 10 -4 
(Savag e~& Sembach]|l996l) . which is also conservative, since the 
gas in the outflow is certainly, on average, at a higher state of 
ionization. The luminosity of the [CII] broad component is in- 
ferred from the flux of this component in the spectrum extracted 
from an aperture of 6" (F[cii](broad) = 6.8 Jy km s -1 ), yielding 
L[cii](broad) = 7.3 + 0.9 x 10' L Q . We obtain a lower limit on the 
outflowing atomic gas mass of 



M outfi (atomic) > 7 x 10 M e 



(2) 



We emphasize that this is a conservative lower limit on the total 
mass of outflowing gas, not only because of the assumptions on the 
physical properties of the outflowing atomic gas, but also because a 
significant fraction of the outflowing gas is likely in the molecular 
form. 



3.3 Extension 

Fig. [2 a) shows the map of the [CII] narrow component integrated 
withinj -300 < v < +400 km s -1 . Note that in this case we 
have subtracted a pseudo-continuum defined by the flux observed 



1 Note that, although Fig.[T]shows the spectrum resampled in channels of 
85 km s - ' for sake of clarity, the fit was performed on the unbinned spec- 
trum to avoid any fitting artifact associated with the binning. 

2 The line core integration limits are asymmetric because the narrow com- 
ponent is slightly skewed towards positive velocities. 




Figure 2. Map of the continuum-subtracted [CII] line narrow compo- 
nent (a), -300 < v < +400 km s~\ and of the [CII] line wings (b), 
-1300 < v < -300 km s -1 and +400 < v < +1300 km s -1 . See text for 
details on the continuum subtraction of the narrow component. The dashed 
circles indicate the extraction apertures of the two spectra shown in Fig.[TJ 
Levels are in steps of 0.64 ly km s -1 beam -1 (i.e. 3tr) in the narrow com- 
ponent map (a) and in steps of 0.36 Jy km s -1 beam -1 (i.e. lcr) in the broad 
wings map (b). The beam of the observation is shown in the bottom-left cor- 
ner. At the redshift of the source 1" corresponds to 5.5 kpc. The colorbars 
are in units of Jy beam -1 , giving the average flux density in each velocity 
integration interval; to obtain fluxes in units of Jy km s -1 beam -1 one has 
to multiply by the width of the integration interval of each map. 



at 400 < |v| < 800 km s -1 (resulting in a level of 5.6 mJy), to mini- 
mize the contribution of the broad component. Fig.|2jb) shows the 
map of the [CII] broad wings, where we have combined the blue 
(-1300 < v < -300 km s -1 ) and red (+400 < v < +1300 km s -1 ) 
wings to improve the signal-to-noise. 

The narrow component (Fig|2^), tracing gas and star forma- 
tion in the host galaxy is resolved at high significance. To illustrate 
this extension more quantitatively, Fig. [3^ shows the amplitude of 
the visibilities as a function of the baseline distance (black sym- 
bols) for the narrow component of the [CII] line. Errorbars show 
the statistical photon noise on the average amplitude in each bin. In 
this diagram an unresolved source has a flat distribution. Such a flat 
distribuion can be ruled out by the uv data. A fit of the uv data with a 
simple circular Gaussian gives a size FWHM=1.5", corresponding 
to 8 kpc. However, as illustrated by the red symbols, showing the 
result of the fit, a single circular Gaussian does not provide a good 
description of the data, indicating that the morphology of the [CII] 
core emission is more complex. We know that the [C II] core has 
a very compact component (size 0.3", Walter et al. 2009) which is 
unresolved at our resolution. This compact component is responsi- 
ble for the flat visibilities at baselines longer than 40 m. The rising 
visibilities at baselines shorter than 40 m indicate the presence of a 
resolved component, with size larger than 2", whose flux is about 
40% of the total narrow line flux. It is likely that this extended com- 
ponent of the narrow [CII] emission was resolved out in previous 
high angular resolution observations. 

The spatial extension of the line wings is more difficult to 
quantify due the lower signal-to-noise. The map of the wings 
(Fig[2p) is characterized by a compact core, detected at 8<x, and 
diffuse emission surrounding the source on scales of a few arcsec, 
but at low significance. The amplitude versus wv-radius diagram, 
shown in Fig.|3p, suggests that the the outflow is indeed marginally 
resolved, as indicated by the decreasing visibilities out to 70 m. The 
visibilities at distances <70m are inconsistent with a flat distribu- 
tion at confidence level higher than 99%. A fit of the uv visibilities 
with a simple circular Gaussian (shown with the red symbols) gives 
a size of 2.9" (16 kpc). If confirmed (with higher resolution and 
deeper observations) this would be the largest outflow ever detected 



© 2002 RAS, MNRAS 000, [Tp?? 



4 Maiolino et al. 



Figure 3. Amplitude of the visibilities as a function of the baseline distance 
(black symbols) for the narrow component of the [CII] line (a) and for its 
broad wings (b). In both cases the distribution of visibilities is inconsistent 
with an unresolved source (which would give a flat distribution). The reds 
symbols show the result of the best fit by using a circular Gaussian model. 
Note that the units (Jy) refer to the average flux in the selected channels 
associated with the line, to obtain the units in terms of line flux one has to 
multiply by 700 km s -1 and by 1800 km s , respectively. 



and it would support models ascribing the capability of cleaning en- 
tire galaxies, out to large radii, to quasar-driven outflows. However, 
the intensity of the last visibility at 90 m suggests that the structure 
of the outflow is more complex and cannot be modelled with a sim- 
ple Gaussian. We have attempted to fit a ring-like morphology to 
the uv data, but the fit does not improve and the inner radius is set 
to a very small value (<0. 1"). 

Finally, we note that the fit also requires an offset of the core 
of the wings of 0.3" (±0. 1") relative to the line core, this is also ob- 
served in the map (Fig.|2p). The offset is small and marginally sig- 
nificant. However, even if real, an offset of the centroid of the out- 
flowing gas is not surprising, since outflows observed i n local and 
high-z galaxies are not aligned with the galaxy centers dWeifi et alj 
l200ll ; IWalter et alj2004 IWilson et alj|2008l : lAlatalo et alfeoill) . 



3.4 Outflow rate, energetics and depletion time 

Given the limited information available, we cannot elaborate on 
sophisticated models. We have assumed a simple scenario where 
the outflow occurs in a spherical volume (or conical or multi- 
conical volume), with radius R ol , t flow> uniformly filled with out- 
flowing clouds. Note that the outflowing clouds can have any den- 
sity distribution, since we are assuming the conservative case of 
n » n cr , which gives a lower limit on the outflowing mass accord- 
ing to Eq. 1. The presence of clouds with density below the critical 
density would result into an even higher outflow mass and outflow 
rate. 

For the geometry discussed above, the volume-averaged den- 
sity of atomic gas in the outflow is (p olll n)v = 3 M oul n/(£2R^ utfl ) 
where CI is the total solid angle subtended by the (multi-)conical 
outflow (CI = An for a spherical outflow). Note that (p ollt fl)v should 
not be confused with the gas density of the individual clouds. In the 
case of a shell (or fragmented shell) geometry the volume-averaged 
density is (p outfl ) v = 3 /^(R^-R^)] * M olltfl /(n R 2 mm dR), 
where Rj n is the inner radius of the shell and where the last approx- 
imation holds in the case that the shell thickness (dR = R ou1 h - Ri„) 
is much smaller than its radius. In these equations CI is the total 
solid angle subtended by the (fragmented) shell. 

In the case of a spherical, or multi-conical, geometry the out- 
flow rate is then given by 



Mo, 



' ^ Routfl (Poutfl)v - 3 V 



M„, 



where v is the outflow velocity. In the case of a (fragmented) shell- 
like geometry 



M„, 



-j R out fl Moutfl M oulfl 

3v — ; — K V ■ 



R 



R 3 



dR 



(4) 



Clearly the shell-like geometry gives an even higher outflow rate 
relative to the spherical/multi-conical geometry. In the following 
we conservatively assume the latter geometry. Note also that the 
outflow rate is independent of CI, both in the (multi-) conical and in 
the (fragmented) shell geometries. 

Since we have a lower limit on M outll , we can derive a conser- 
vative lower limit on the outflow rate. In particular, by taking R oul n 
= 8 kpc and v = 1300 km s (conservatively assuming that the 
maximum velocity observed in the wings is the de-projected veloc- 
ity of the outflow, which is a lower limit in the case of a conical 
outflow not intercepting the plane of the sky), we obtain 



M™ 



> 3500 M yr 



(5) 



We note that if the outflow size has been overestimated, due to 
the low signal-to-noise in the wings, as discussed above, and the 
bulk of the flux is actually unresolved, then the lower limit on the 
outflow rate in Eq.[5]is even more conservative, since the outflow 
rate scales as 1 /Routfl- As discussed above, a (fragmented) shell- 
like geometry would given a even higher outflow rate. Concerning 
the velocity field, the angular resolution of our data does not al- 
low us to map the kinematics of the outflow. However, we note 
that, in the nearby quasar Mrk231, transitions tracing the outflow 
on differ ent scales have the same profile and the same maximum 
velocity dAalto et ani2012l ; lcicone et alj|2012l) . supporting the sce- 
nario where quasar driven outflows have velocity roughly constant 
throughout the wind. The same kinematics is likely to apply to 
J 1148 as well. However, new higher angular resolution observa- 
tions will allow us to better constrain the kinematics. 

The inferred lower limit on the outflow rate is similar to the 
star formation rate in the host galaxy (~ 3000 M yr -1 ). 

The kinetic power associated with the outflow is given by 



0.5 v 2 M outfl > 1.9 x 10 45 ergs- 



(6) 



(3) 



which is about 0.6% of the bolometric luminosity of the quasar nu- 
cleus. On the one hand, this is well below the maximu m efficiency 
expec ted for quasars to drive outflows (Px/Lboi ~ 5%. lLapi et alj 
2005), meaning that the quasar radiation pressure is probably fully 
capable of driving the outflow observed in Jl 148+5251. On the 
other hand, the lower limit on the outflow kinetic power is barely 
consistent with the kinetic power that can be achieved by a star- 
burst driven wind (P K (SB) = 7 x 10 4 ' (SFR/Mp y r 1 ) erg s -1 « 
2 x 10 45 erg s -1 for Jl 148+5 152. IVeilleux et alj2005l) , therefore we 
cannot exclude a starburst contribution to the wind based solely on 
energetic arguments and by using simple scaling equations. How- 
ever, as discussed above the very high outflow velocities are dif- 
ficult to explain by starburst driven winds models, hence strongly 
favoring quasar radiation pressure as the main driving mechani sm. 

We also mention that a parallel paper (Valia nte et al.l2012h . by 
exploiting existing detailed cosmological models specifically suited 
for Jl 148 constrained by all observable quantities, do predict a 
quasar-driven wind fully consistent with the outflow rate measured 
by us, while the outflow contribution by the starburst is orders of 
magnitude lower. 

The molecular gas content in the host gala xy of Jl 148+5251 , 
as inferred by CO ob servations, is 2 x 10 10 M Q (Walt er et alj2003l: 
Bertol di et al . 2003b). At the observed (minimum) outflow rate, the 
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quasar host galaxy will be cleaned of its gas content, hence quench- 
ing star formation, within less than 6 x 10 6 yr. 



4 CONCLUSIONS 

By using the IRAM PdBI, we have discovered an exceptional out- 
flow in the host galaxy of the quasar Jl 148+525 1 at z=6.4, through 
the detection of prominent broad wings of the [CII]158yum line, 
with velocities up to about +1300 km s -1 . This is the most dis- 
tant massive outflow ever detected and it is likely tracing the long 
sought quasar-feedback already in place in the early Universe, close 
to the re-ionization epoch. Both the high outflow velocity and the 
kinetic power (Pk > 1.9 x 10 45 erg s -1 ) favor quasar radiation pres- 
sure as the main driving mechanism. A more detailed cosmological 
model of Jl 148 dValiante et al]|2012h confirms that in this object 
a quasar-driven outflow is expected with an outflow rate consitent 
with what observed by us, while the contribution from the starburst 
wind is orders of magnitude lower. 

We marginally resolve the outflow, with a size of about 16 kpc, 
which would make it the largest outflow ever detected. This finding 
supports models ascribing the capability of cleaning entire galax- 
ies, out to large radii, to quasar-driven outflows. However, higher 
angular resolution and deeper observations are required to confirm 
the outflow extension. 

We infer a conservative lower limit on the total outflow rate 
of Moutfl > 3500 M yr , which would be the highest outflow rate 
ever detected. This lower limit on the outflow rate is higher than 
the SFR in the quasar host (SFR » 3000 M Q yr" 1 ). At the observed 
outflow rate the gas content in the host galaxy will be cleaned, and 
therefore star formation will be quenched, in less than 6 x 10 s yrs. 
Such a fast and efficient quenching mechanism, already at work 
at z>6, is what is required by models to explain the properties of 
massive, old and passive galaxies observed in the local Universe 
and at z ~ 2. 
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